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We consider the separability of various joint states of D-dimensional quantum sys- 
tems, which we call "qudits." We derive two main results: (i) the separability condition 
for a two-qudit state that is a mixture of the maximally mixed state and a maximally 
entangled state; (ii) lower and upper bounds on the size of the neighborhood of sepa- 
^ ■ rable states surrounding the maximally mixed state for N qudits. 
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1 Introduction 

^ ' One of the distinguishing features of quantum mechanics, not found in classical physics, is 
■ the possibility of entanglement between subsystems. It lies at the core of many applications 
in the emerging field of quantum information science such as quantum teleportation 
and quantum error correction |P, Entanglement is a distinctly quantum- mechanical 
correlation between subsystems, which cannot be created by actions on each subsystem 
separately; moreover, correlations between subsystem measurements on an entangled com- 
posite system cannot be explained in terms of correlations between local classical properties 
inherent in the subsystems. Thus one often says that an entangled composite system is 
nonseparable. Formally, the state of a composite system, pure or mixed, is separable if the 
state has an ensemble decomposition in terms of product states. A separable state has no 
quantum entanglement, and a nonseparable state is entangled. Though the nonclassical na- 
ture of quantum entanglement has been recognized for many years |^ ^ , only recently has 
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considerable attention been focused on trying to understand and characterize its properties 
precisely. 

This paper focuses on the question of whether various joint quantum states of D-dimen- 
sional quantum systems are entangled. For convenience, we call a D-dimensional quantum 
system a "qudit," by analogy with the name "qubit" for D = 2 and "qutrit" for D = 3. We 
now have a general method for quantifying the degree of entanglement of a pair of qubits [0, 
and we have a criterion, the partial transposition condition of Peres ||^, which determines 
whether a general state of two qubits is entangled and whether a general state of a qubit 
and a qutrit is entangled 0. The partial-transposition condition fails, however, to provide a 
criterion for entanglement in other cases, where the constituents have higher Hilbert-space 
dimensions |]TD|, |TT| or where there are more than two constituents. Indeed, at present there 



is no general criterion for determining whether the joint state of qudits is entangled, nor 
is there any general way to quantify the degree of entanglement if such a state is known to 
be entangled. 

In Sect, ^we review an operator representation of qudit states, which is applied in Sect. ^ 
where we consider states of two qudits that are a mixture of the maximally mixed state 
and a maximally entangled state. We show that such states are separable if and only if the 
probability for the maximally entangled state in the mixture does not exceed 1/{1 + D). This 
result was obtained by Horodecki and Horodecki [1^, and a more general result, of which 
this is a special case, was obtained by Vidal and Tarrach [|T^. In Sect. ^ we consider the 
separability of mixed states of qudits near the maximally mixed state. We find both lower 
and upper bounds on the size of the neighborhood of separable states around the maximally 
mixed state. Our results generalize and extend the results obtained by Braunstein et al. for 
qubits and by Caves and Milburn for qutrits |jl5|. Before tackling the upper and lower 



bounds, we present, in Sect. |]T, various mathematical results which are used to obtain the 
lower bound, but which might prove useful in other contexts as well. 



2 Operator representation of qudit states 

In this section we review an operator representation of qudit states, analogous to the Pauli, 
or Bloch-sphere, representation for qubits. We begin with the set of Hermitian generators of 
SU(D); the generators, denoted by Aj, are labeled by a Roman index taken from the middle 
of the alphabet, which takes on values j = 1, . . . , — 1. We represent the generators in an 
orthonormal basis |a), labeled by a Roman letter taken from the beginning of the alphabet, 
which takes on values a = 1, . . . , Z). With these conventions the generators are given by 

j = l,. 

A, = r.^ ^ / (Y,\b){b\-{a-l)\a){a\ \ , 2<a<D, (1) 
aia — 1) \b=i / 



J = /},..., (D + 2)(D-l)/2: 

A, = ri:) = -^(|a)(6| + |6)(a|) , l<a<h<D, (2) 
3 = D{D + 2)/2,...,D^-l: 

A,=rl^) = ^(|a)(6|-|6)(a|) , l<a<h<D. (3) 
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In Eqs. (U) and @, the Roman index j stands for the pair of Roman indices, ab, whereas in 
Eq. (|ip, it stands for a single Roman index a. The generators are traceless and satisfy 

XjXk = '^^jk + djkiXi + ifjkiXi . (4) 

Here and wherever it is convenient throughout this paper, we use the summation convention 
to indicate a sum on repeated indices. The coefficients fjki, the structure constants of 
the Lie group SU(Z)), are given by the commutators of the generators and are completely 
antisymmetric in the three indices. The coefficients djki are given by the anti-commutators 
of the generators and are completely symmetric. 

By supplementing the — 1 generators with the operator 

Ao = , (5) 



where / is the unit operator, we obtain a Hermitian operator basis for the space of linear 
operators in the qudit Hilbert space. This is an orthonormal basis, satisfying 

tr(A„A/3) = 5af3 ■ (6) 

Here the Greek indices take on the values 0, . . . , — 1; throughout this paper, Greek indices 
take on or more values. Using this ort honor mality relation, we can invert Eqs. 
to give 

= ^+ / / , |-(a-i)r,+ E r^l , (7) 

^ \ aia — 1 




\a){b\ = —(rit^ + zri^)), l<a<b<D, (8) 

\b){a\ = -^(rW-^ri^)), l<a<b<D. (9) 
Any qudit density operator can be expanded uniquely as 

p = ^CaXa , (10) 

where the (real) expansion coefficients are given by 

c„ = DtT{pX^) . (11) 
Normalization implies that Cq = V^, so the density operator takes the form 

P = ^iI + c,X,) = ^iI + c-X). (12) 

Here c = CjCj can be regarded as a vector in a (D^ — l)-dimensional real vector space, spanned 

by the orthonormal basis ej, and A = XjCj is an operator-valued vector. If p = is a 
pure qudit state, then tr(p^) = 1, from which it follows that 

\c\'^ = c-c = D{D-1) . (13) 



We could represent a pure state by a unit vector n = c/ y D{D — 1) on the unit sphere in 
— 1 dimensions, but in contrast to the situation with the Bloch sphere {D = 2), most 
vectors on this unit sphere do not represent a pure state or, indeed, any state at all. 
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3 Mixtures of maximally mixed and 
maximally entangled states 

In this section we deal with two qudits, labeled A and B. We consider a class of two-qudit 
states, specifically mixtures of the maximally mixed state, = I^I/D"^, with a maximally 
entangled state, which we can choose to be 

I^) = 4hEI«)®I«)- (14) 

V-t^ a=l 



Such mixtures have the form 

= (l-e)Mz52+e|^)(^| , (15) 

where < e < 1. 

In analogy to Eq. (p!OD, any state p of two qudits can be expanded uniquely as 

P = -^CapK ® A/3 , (16) 

where the expansion coefficients are given by 

= Dhi{pK ® , (17) 

with Coo = -D determined by normalization. Using Eq. (|l^ or Eqs. (0)-(0), we can find the 
operator expansion for the maximally entangled state (p!^ : 

1^) (^1 = 7^ M ® ^ + ^ E Ta ® + D 5: (rw ® - rl^ ® rl, ] , (is) 

^ \ a a<b J 

from which we can read off the expansion coefficients for the state pe of Eq. (|T3|): 

(19) 

j ^ 

j = k = l,...,{D + 2){D-l)/2. (20) 
j = k = D{D + 2)/2,...,D^-l. 

A state of the two qudits is separable if it can be written as an ensemble of product 
states. In this section we show that the mixed state state (|15|) is separable if and only if 

e < . (21) 

- 1 + D ^ ' 

Our method is to prove the necessity of the condition ( PT| ) by considering the restrictions 
that separability places on the correlation coefficients (|20| ) and then to construct an explicit 
product ensemble when e < 1/(1 + Z)). Vidal and Tarrach |jl3| found the separability bound- 
ary for a mixture of with any pure state by using the partial transpose condition |Q 
to show that any state with e outside the boundary is nonseparable and by constructing an 
explicit product ensemble for states with e within the separability boundary. Horodecki and 
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Horodecki [0 found the separability boundary for the state ( p!5D using other techniques. 
The reason for presenting in this section a more hmited result than that of Vidal and Tar- 
rach is, first, that our proof of necessity has a nice physical interpretation in terms of the 
correlation coefficients ( [2D| ) and, second, that the product ensemble we use is different from 
the one used by Vidal and Tarrach. 

The product pure states for two qudits, |'?/'a)('^a|®|'^b)(V'b|, constitute an overcomplete 
operator basis. Thus we can expand any two-qudit density operator in terms of them. 



p= dVAdVBw{lpA,1pB) \i'A){lpA\®\i'B){lp 



B\ 



(22) 



Here the integral for each system runs over all of projective Hilbert space, i.e., the space 
of Hilbert-space rays, and the volume elements (IVa and (IVb are the unitarily invariant 
integration measures on projective Hilbert space. Because of overcompleteness of the pure- 
state projectors, the expansion function w{ipA,'ipB) is not unique. Notice that the expansion 
coefficients Cap of Eq. ([T7|) can be written as integrals over the expansion function. 



Ca/3 = / dVAdVBw{'4)A,i'B){cA)a{cB)p 



(23) 



where (ca)^ = D {ip A\'^a\''P a) and (cb)^ = D{ipB\'^a\4'B) are the expansion coefficients for the 
pure states lip) a and \'iP)b, satisfying ca- ca = D{D — 1) = cb ■ cb- 

A two-qudit state is separable if and only if there exists an expansion function w{iIja, '4'b) 
that is everywhere nonnegative. In this case w{iI)Ai'4^b) can be thought of as a normalized 
classical probability distribution for the pure states ipA and ipB-, and the integral for in 
Eq. (^) can be interpreted as a classical expectation value of the product of the random 
variables {cA)a and (c^)/?, i.e.. 



Ca/3 = E [(cA)a(cB)/3] . 

If the state is separable, we have from Eq. (^) that for each value of j, 

1 



Dt 



E[icA),icB) 



< 



E[{ca)] + E[{cb) 



Adding over the — 1 value of j gives 



D{D^ - l)e < Ue[ca ■ ca] + E[cb ■ cb]) = D{D - 1) . 



(24) 



(25) 



(26) 



We conclude that if pe is separable, then e < + D). 

To prove the converse, we construct an explicit product ensemble for the state p^ with 
e = 1/(1 + D). We define a vector z = {zi, . . . , zd) whose components Za take on the values 
±1 and ±2, so that 



Zj Zj Zj 

Associated with each vector z is a pure state 

1 



(27) 



D 



■ 



{21 



a=l 
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There are 4^ vectors and thus that many states although only 4^~^ of these states are 
distinct in that they differ by more than a global phase. Now we define a product state for 
the two-qudit system: 

PS = I ($,-1 ® I ^g, ) I . (29) 

The ensemble consisting of all 4^ of these states, each contributing with the same probability, 
produces the density operator 



z a,b,c,d \ z / 



(30) 



Since 



it follows that 



ZaZlzlZd = 4:^{SabScd + SacSbd - ^ah^cd^ac) , (31) 



^Ep^^ = ^ + - E i«)H ® • (32) 

z a=l 

Multiplying hj D/{D + 1) and rearranging yields 

The left-hand side of Eq. (^) is the state (|15D with e = 1/(1 + D), and the right-hand side 
is an explicit product ensemble for the state. This concludes the proof that is separable 
if and only if e < l/{l + D). 



4 Separability of states near the 
maximally mixed state 

This section deals with A^-qudit states of the form 

= (1 - e)M£,iv + epi 



(34) 



where Mj^n = J® ■ ■ - ^I/D^ is the maximally mixed state for qudits and pi is any A^-qudit 
density operator. We establish lower and upper bounds on the size of the neighborhood of 
separable states surrounding the maximally mixed state. In particular, we show, first, that 
for 



e < 



1 + D2iV-l ' 



all states of the form (|3^) are separable and, second, that for 

1 



e > 



1 + D^-i ' 



(35) 



(36) 



there are states of the form (0) that are not separable (i.e., they are entangled). These 
results generalize and extend the work of Braunstein et al. for qubits [Q and of Caves and 



Milburn for qutrits |15 . 
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4.1 Mathematical preliminaries 

Before turning to the lower and upper bounds, it is useful to develop some mathematical 
apparatus that will be used in deriving the bounds. 

4.1.1 Superoperator formalism 

We begin by reviewing a formalism for handling superoperators, introduced by Caves |]IB| 



and used by Schack and Caves to generate product ensembles for separable iV-qubit 
states. 

The space of linear operators acting on a D-dimensional complex vector space is a D^- 
dimensional complex vector space. In this space we introduce operator "kets" |A) = A and 
"bras" {A\ = A\ distinguished from vector kets and bras by the use of smooth brackets. 
The natural operator inner product can be written as = tT[A^B). An orthonormal 

basis I a) induces an orthonormal operator basis, 

\c) {a\ = Tea = Ta , (37) 

where the Greek index a is an abbreviation for the pair of Roman indices, ca. Not all 
orthonormal operator bases are of this outer-product form. 

The space of superoperators, i.e., linear maps on operators, is a D^-dimensional complex 
vector space. Any superoperator S is specified by its "matrix elements" 

S,a,db= {c\Si\a){b\)\d) , (38) 

for the superoperator can be written in terms of its matrix elements as 

<S = '^ca,db\c){a\ \b){d\ = J2 <Sca,dbna © T^b = XI '^"/3 1 (^/3 1 " (^9) 

c,a,d,b c,a,d,b a,l3 

The tensor product here is an ordinary operator tensor product, but we use the symbol 
to distinguish it from a tensor product between objects associated with different systems, 
which is denoted by ®. In the final form of Eq. (^), the tensor product is written as an 
operator outer product, with a = ca and P = dh. 

The ordinary action of S on an operator A, used to generate the matrix elements, is 
obtained by dropping an operator A into the center of the representation of iS, in place of 
the tensor-product sign, 

S{A) = Y,S^pt^AtI . (40) 

There is clearly another way that S can act on A, the left-right action, 

5|A) = X5.^|r.)(r^|A), (41) 

in terms of which the matrix elements are 

So.p = (r,| S\Tp) = (r,,| S\ub) = (c| S(\a){h\)\d) . (42) 

This expression provides the fundamental connection between the two actions of a superop- 
erator. We can define an operation, called sharp, that exchanges the ordinary and left-right 
actions: 

S*{A) = S\A) . (43) 
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noS= n^sS^pT^T^Qrlrl . (48) 



Equation ([42|) implies that 

Sfa,,, = (c| S*{\a){b\)\d) = [rJSK,) = 5,,,,, (44) 

or, equivalently, that 

S*= Sca,db\c){d\Q\b){a\ . (45) 

c,a,d,b 

With respect to the left-right action, a superoperator works just like an operator. Mul- 
tiplication of superoperators TZ and S is given by 

nS = Y '^ay<S^t3\ra)irp\ , (46) 

and the adjoint is defined by 

{A\S^\B) = iB\S\Ar ^ S^ = T.^;Mir,\. (47) 

a,l3 

With respect to the ordinary action, superoperator multiplication, denoted as a composition 
71 o S, is given by 

The adjoint with respect to the ordinary action, denoted by , is defined by 

tT{[S^B)U)=tT{B^S{A)) ^ 5x=E5:,rt0r,. (49) 

The identity superoperator with respect to the left-right action can be written as 

I = El^°)(^"l =Elc)(a|0|«)(c| • (50) 

a c,a 

When sharped, I becomes the identity superoperator with respect to the ordinary action, 
denoted by I: 

I* = Y,\c){c\Q\a){a\ = IQI =1 . (51) 

c,a 

The final ingredient we need is the superoperator trace relative to the left-right action, 
defined by 

Tr(5) = j:i^Mra) = E(c| '5(|a)(a|)|c) = tr(5(/)) . (52) 

a c,a 

Notice that I(/) = DI and J(/) = /, which give Tr(I) = and Tr(X) = D. 

Now suppose the operators \Na) constitute a complete or overcomplete operator basis; 
i.e., let the operator kets \Na) span the vector space of operators. It follows that the super- 
operator Q defined by 

G = J2\Na)iN^\ = Q^ (53) 

a 

is invertible with respect to the left-right action. The operators 

\Qa)=G-'\N^) (54) 



form a dual basis, which gives rise to the following expressions for the identity superoperator: 

I = \Qa){Na\ = E \Na){Qa\ • (55) 
a a 

An arbitrary operator A can be expanded in terms of the original basis or the dual basis: 

^ = T.\No.){Qa\A) = Y.NMQiA) , (56) 

Q a 

A = J2\Qc.)iNa\A) = J2QMNiA). (57) 

a a 

These expansions are unique if and only if the operators \Na) are linearly independent. Later 
in this section we apply expansions of this sort to density operators. 



4.1.2 Pure states and their dual basis 

The set of all pure-state projectors in a D-dimensional Hilbert space, 

P^ = \i;){i;\, (58) 

forms an overcomplete operator basis. To develop operator expansions in terms of the 
pure-state projectors, we follow the discussion in the preceding subsection and consider the 
superoperator 

g = JdV\P^){P^\= JdV\ij){ij\Q\^){^\, (59) 

where dV is the unitarily invariant integration measure on projective Hilbert space. 

The only Hilbert-space integrals we need to calculate explicitly are those for which the 
integrand is a function only of an angle 6 defined by cos^ = KelV^)!, where |e) is some 
particular unit vector (pure state). The angle 6, which runs over the range < ^ < 7r/2, can 
be thought of as a "polar angle" relative to the "polar axis" defined by |e). For integrals of 



this sort, a convenient form of the integration measure is 

dV = (sin 6)^^-^ cos 6 dO dSio-z , (60) 
where dS2D-3 is the standard integration measure on a {2D — 3)-dimensional unit sphere. 



Thus the total volume of D-dimensional projective Hilbert space is [18 



V = S,n^s ['^ de (sin ef^~' cos 6 dO = = , (61) 

where S2D-3 = 27t^~^/{D — 2)! is the volume of a {2D — 3)-dimensional unit sphere. 

To use the expansions (|56| ) and (^Tf), we need the dual basis \Q^), and for that purpose, 
we need to invert Q. Since Q is Hermitian relative to the left-right action, we can invert it 
by diagonalizing it with respect to the left-right action. Given an orthonormal basis |a), we 
can write Q as in Eq. (p9D, 



= Gca,db\c){a\&\b){d\ = Gca,db\Tca){rdb\ , (62) 

c,a,d,b c,a,d,b 
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where the matrix elements are given by Eq. (|38D : 

Qca,db = {c\g{\a)md) = I dV{c\^){iP\a){h\^){m ■ (63) 



The unitary invariance of the integration measure places stringent constraints on the 
matrix elements (]63|). Since the integral in Eq. ( |63D remains unchanged under a change in 
the sign of the amplitude {a\il)) corresponding to a particular basis vector |a), the matrix 
elements vanish except when {})a = h^c = doia = c^h = doi (ii) a = h = c = d. 
Furthermore, unitary invariance implies that for each of these cases, all the matrix elements 
have the same value. Gathering these conclusions together, we have 

{a, a = b^c = doTa = c^b = d, 
7, a = b = c = d, (64) 
, otherwise. 

We get a relation between a and 7 by noting that 

D 

D{D -l)a + D^= J2 Gca,ca = V , (65) 

c,a=l 



where the second equality follows from doing the sum within the integral in Eq. (|63|). We 
need one more relation, which we get by evaluating explicitly the integral for 7: 



J dV\{a\^)\' = S2D-3 I de{Mf'>-%cosef = ^^^—^^2K . (66) 



7 

It follows that a = K. As a result, we have 



g = K(2Y^\Taa){Taa\ + E {Wb){rab\ + \raa){n,\) \ = K{l + I) . 
\ a a,b / 



(67) 



This result gives us immediately that |T8 



ldV\^){^P\ = g{I) = ^I . (68) 

The operators Aq, introduced in Sect. ^ constitute a complete, orthonormal operator basis, 
so we can write I as 

I = ElA«)(A«l = ^ + ^> (69) 

where 

r = E|A,)(A,| (70) 
i 

is the superoperator that projects onto the subspace of traceless operators. Plugging Eq. ( |5U[ ) 
into Eq. (^) gives the diagonal form of Q: 

6^ = i^(^(/^ + i)^ + r^ . (71) 
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Orthonormal eigenoperators of Q are Aq = I /\fD, with eigenvalue K{D + 1) = V/-D and the 
traceless operators Aj, which are degenerate with eigenvalue K = V / D{D + 1). 

We are now prepared to write the inverse of Q with respect to the left-right action as 

Thus the dual operators of Eq. (^ID are given by 

m=Q-'\P^) = ]^(\P^)-^^=^{{D + l)P^-l) . (73) 

4.1.3 Alternative diagonalization of Q 

In this subsection we rederive Eq. ( fTTj ) using the special properties of the superoperator Q. 
These properties are evident from the integral form of Q in Eq. (|59D . 

• The superoperator Q is Hermitian relative to the left-right action, which implies that 
it has a complete, orthonormal set of eigenoperators r/^, a = 1,...,D^, with real 
eigenvalues qa- 

a a 



• The superoperator Q is Hermitian relative to the ordinary action, 

^ = = E ic^vi qv.=y: qMMi , (75) 

a a 

which implies that if rja is an eigenoperator of Q, then rjl^ is also an eigenoperator 
with the same eigenvalue. This means that we can choose all the eigenoperators to be 
Hermitian. 

• The superoperator Q is unitarily invariant, i.e., 

g = UQU^ ogou^ QU = Y.iaUvaU^ & u^lu^ , (76) 

a 

for any unitary operator U, which implies that if rja is an eigenoperator of Q, then 
UrjoiU'^ is also an eigenoperator with the same eigenvalue. 

The upshot of these three properties is that the eigensubspaces of Q are invariant under 
Hermitian conjugation and under all unitary transformations. It is not hard to show that 
the only such operator subspaces are the subspace of traceless operators and its orthocom- 
plement, the one- dimensional subspace spanned by the unit operator. The result is that Q 
must have the form 

a = A-(^^ + T)=4 + ^i^l). (77) 

where K is the eigenvalue of any traceless operator and K^i is the eigenvalue of Aq = I /\fD. 
Now we use the final property to evaluate /i. 
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• The superoperator Q is invariant under exchange of the two kinds of action: 

g = g# = K{l+^^i^ . (78) 

This imphes that yU = D + 1, thus bringing Q into the form (p7|), but with K not yet 
determined. 

We find the value of K by evaluating the superoperator trace, first using Eq. (^), 

Tr(^) = tr(^(/)) = V , (79) 

and then using Eq. (|67D, 

Tr(^) = k{Ti{1) + Tr(J)) = KD{D + 1) . (80) 
This gives K = V / D{D + 1), in agreement with Eq. (|66D . 

4.2 Separability bounds 

We turn now to demonstrating the lower and upper bounds, Eqs. ( p5D and (|36D, on the size 
of the neighborhood of separable states surrounding the maximally state. 

To establish the lower bound, we use the results of Sect. ^T!] to formulate operator ex- 
pansions in terms of product pure states. For a single qudit, any density operator can be 
expanded as 

p = [ dV\P^){Q4p) = I dVw,{^)P^ , (81) 



where 

w.iij) = tr(pg^) = ^{iD + imp\^) - l) (82) 

is a quasi-probability distribution, normalized to unity, but possibly having negative values. 
The analogous product representation for an A^-qudit density operator is 

p = J dVi--- dVN Wpiijji, . . . , i/jn) P^i ® ■ ■ ■ ® Pi>N ) (83) 

where 

Wp{ipi, iPn) = tr(p(5^^ (g) ■ ■ ■ ® Q^j^) . (84) 
The A^-qudit quasi-distribution obeys the bound 

/ , , \ ^ / smallest eigenvalue of \ jj'^n-i 

This follows from the fact that has a nondegenerate eigenvalue, /V, and a.{D — l)-fold 
degenerate eigenvalue, —D/V. Thus the most negative eigenvalue of the product operator 
Q^, ® ■ ■ ■ ® Qv'iv is (-D/V)(DVV)^-i = -D^A^-i/V^. 

We can use the lower bound ( P5| ) to place a similar lower bound on the quasi-distribution 
for the mixed state p^ of Eq. (|3^). Since the quasi-distribution for the maximally mixed 
state, M^N, is the uniform distribution 1/V^, we have 

1-e 1 - e(l + Ll^iv-n 
Wp,(^i,...,V'7v) = ^j^ + ewpi > y]^ • (86) 
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We conclude that if e < 1/(1 + D"^^^^), then Wp^ is nonnegative and the qudit state is 
separable. This estabhshes the lower bound ( P^D on the size of the neighborhood of separable 
states surrounding the maximally mixed state. 

The upper bound (^) on the size of the separable neighborhood can be established with 
the help of an exact separability condition for a particular A^-qubit state, obtained by Diir, 
Cirac, and Tarrach |]l9l and also by Pittenger and Rubin pO|. We consider the iV-qudit 



state, 

= (1 - e)MDN + e|^cat)(^cat| , (87) 

where 

1 ^ 

l^cat) = l«) ' (88) 

a=l 

is an A^-qudit "cat state." We call the mixed state (0) an e-cat state. 

Now project each qudit onto the two-dimensional (qubit) subspace spanned by |1) and 
|2). The local projection operator on each qudit is 11 = |1)(1| + |2)(2|, and the normalized 
iV-qubit state after projection is 

P'e = ^.(nJv) = (1 - + e1$cat)($cat| , (89) 

where 

|$eat) = -1={\1) ® . . . ® |1) + |2) ® . . . ® |2)) (90) 



V2 

is the cat state for qubits and 



{2/D)^{l-e) + 2e/D 



Diir, Cirac, and Tarrach and also Pittenger and Rubin |^ have shown that the A^-qubit 
e-cat state (^) is nonseparable (entangled) if and only if e' > 1/(1 + 2^^^), a condition 
equivalent to e > 1/(1 -|- D^~^). Since local projections on each qudit cannot create entan- 
glement, we can conclude that the A^-qudit e-cat state ( ^T]) is nonseparable under the same 
condition. This establishes the upper bound (|36D on the size of the separable neighborhood 
around the maximally mixed state. 



Pittenger and Rubin |^ have recently extended the result of Diir, Cirac, and Tarrach 
for the A^-qubit e-cat state. They have shown directly that the A^-qudit e-cat state (^7^ 
is nonseparable if e > 1/(1 + D^~^), and they have also shown that the same condition is 
a necessary and sufficient condition for entanglement when D is prime. Their argument is 
akin to the correlation-coefficient argument we give in Sect. & 
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